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ABSTRACT 


Each  member  of  the  alkali  metal  series  has  been  investigated  to  determine  the 
best  system  for  a  closed  cycle  MHD  generator  utilizing  magnetically  induced  non¬ 
equilibrium  ionization.  MHD  generator  performances  have  been  calculated  for 
one-dimensional  steady  state  flows  with  constant  Mach  number  in  which  the  dimer 
concentration  is  neglected.  The  calculations  are  for  a  generator  with  segmented 
electrodes. 

Potassium  appears  to  be  the  best  choice  among  the  alkali  metal  series  for  a 
1  megawatt  generator  operating  at  total  temperatures  near  1600®K  and  total  pressures 
around  20  psia.  This  choice  is  based  upon  considerations  of  the  calculated  generator 
performance,  the  vapor  pressure,  and  the  equilibrium  concentration  of  dimer  in  the 
vapor. 

The  calculations  show  that  the  performance  of  an  MHD  generator  utilizing  a 
non-equilibrium  condition  of  the  electrons  is  primarily  dominated  by  the  elastic 
electron-neutral  collision  cross  section.  If  the  elastic  electron  collision  cross  section 
in  potassium  is  4x10  cm^,  then  magnetic  fields  of  130,000  gauss  will  be  required 
to  obtain  the  non-equilibrium  conditions  necessary  to  obtain  power  densities  near 
1  kw/cm'^  when  the  Mach  number  is  around  2.  However,  if  the  cross  section  in 
potassium  is  near  1x10  ^  cm^,  (as  recently  reported  in  the  literature)  then  magnetic 
fields  of  12,000  gauss  will  produce  power  densities  near  1  kw/cm^  at  a  Mach  of  2. 
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Introduction 


In  order  to  operate  a  closed  loop  MHD  generator  with  high  power  density  at 
temperatures  compatible  with  existing  ’’high  temperature  materials"^  a  non¬ 
equilibrium  condition  of  the  electrons  must  be  created  in  the  plasma.  Based  upon 
the  theory  of  magnetically  induced  non-equilibrium  ionization  (Ref.  5),  calcula¬ 
tions  have  been  performed  for  each  member  of  the  alkali  metal  series.  The 
calculations  are  based  upon  a  flow  which  is  one-dimensional,  steady  state  and 
in  which  the  dimer  concentration  is  negligible. 

These  calculations  together  with  considerations  of  the  vapor  pressure  and 
equilibrium  dimer  concentration  indicate  that  potassium  is  the  best  choice  of 
working  fluids  among  the  alkali  metal  series. 
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Equations! 

For  the  case  of  segmented  electrodes,  the  electron  temperature  is  related  to 
the  total  gas  temperature  by  the  following  equation,  (Ref.  5); 


i.-  >♦'3  'f  o-K)'' 


where 


Tg  is  the  electron  temperature, 

Tq  is  the  total  gas  temperature. 

y  is  the  ratio  of  heat  capacities,  (cp  /  c^). 

CO  is  the  electron  cyclotron  frequency, 

is  the  average  time  between  electron-non~electron 
col  lisions . 


is  Mach  number. 

^  is  the  loss  factor  which  is  close  to  unity  for  vapors  with 
monatomic  particles. 

K  is  the  loading  factor  which  is  the  ratio  of  the  load  voltage 
to  the  open-circuit  voltage. 


The  electron  temperatures  are  calculated  for  the  case  where  the  Mach  number  is 
constant,  ani  where  if*  =  1  and  K  =  0. 
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The  electron  density  is  determined  from  the  Sc^ia  equation  based  upon  the 


electron  temperature: 

■We’  _  -  T::^) 


(2) 


where 

n  is  the  electron  density, 
n^  is  the  neutral  particle  density, 
is  the  electron  mass, 
k  is  the  Boltzmann  constant, 
h  is  the  Plank  constant. 

Eq  is  the  first  ionization  potential  of  the  vapor. 

Tg  is  the  electron  temperature. 

The  neutral  particle  density  is  determined  from  the  ideal  gas  law: 


where 

n  ^  is  the  neutral  particle  density, 
P'1  is  the  static  pressure. 

T|  is  the  static  temperature, 
k  is  the  Boltzmann  constant. 
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The  static  pressure  Is  calculated  from  the  isentropic  flow  equation: 


where 

is  the  static  pressure, 
is  the  total  pressure. 

y  is  the  ratio  of  heat  capacities,  (cp  /c^). 
M-j  is  the  Mach  number. 


The  static  temperature  is  calculated  from  the  isentropic  flow  equation: 


T-  T 


(5) 


where 


y 


is  the  static  gas  temperature. 

is  the  total  gas  temperature. 

is  the  static  gas  pressure. 

is  the  total  gas  pressure. 

is  the  ratio  of  heat  capacities,  (Cp/ c^)* 
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The  conductivity  due  to  the  ions  is  related  to  the  electron  temperature,  the 
Debye  shielding  length,  and  the  average  impact  parameter  by  the  following 
equation,  (Ref.  10): 


I 


l.st  KIO^ 

A 


(6) 


where 

is  the  conductivity  due  to  the  ions. 

T  is  the  electron  temperature, 

^  is  the  ratio  of  the  Debye  shielding  length  to  the  average 
impact  parameter. 


The  ratio  of  the  Debye  shielding  length  to  the  average  impact  parameter 
is  related  to  the  electron  density  and  electron  temperature  by  the  following 
equation,  (Ref.  10): 


A 


3 


\  Trrtg 


) 


(7) 


where 


Z  is  the  number  of  electronic  charges  associated 
with  the  ion , 

is  the  number  of  electronic  charges  associated 
with  the  electron . 

e  is  the  electron  charge. 

k  is  the  Boltzmann  constant. 

Tg  is  the  electron  temperature. 

is  the  electron  density. 
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Since  the  case  of  singly  ionized  particles  is  considered,  Z  =  ==  1 . 

The  conductivity  due  to  the  neutral  particles  is  determined  from  the  scalar 
conductivity  equation  in  which  the  electrons  are  assumed  to  have  a  Maxwellian 
distribution: 


where 

e  is  the  electron  charge. 

m  is  the  electron  mass- 
e 

k  is  the  Boltzmann  constant, 
n^  is  the  electron  density. 

Tg  is  the  electron  temperature, 
n^  is  the  neutral  particle  density, 
n^  is  the  electron  density. 

Qenis  electron-neutral  elastic  collision  cross  section. 


The  calculations  are  performed  with  a  constant  electron-collision  cross 
section,  thus  neglecting  the  Ramsauer  effect. 

The  plasma  electrical  conductivity  is  determined  from  the  relation: 


-.1 


(9) 
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The  gas  velocity  is  determined  from  the  sonic  velocity  equation: 


where 


y 

R 


m 

n 


is  the  gas  velocity . 

is  the  ratio  of  heat  capacities,  (c  /c  ). 

p  V 

is  the  universal  gas  constant, 
is  the  static  temperature, 
is  the  atomic  weight  of  the  vapor, 
is  the  Mach  number. 


The  relation  between  the  magnetic  field,  the  electron  density,  the  electron 
cyclotron  frequency,  the  average  time  between  electron-non-electron  collisions, 
and  the  plasma  electrical  conductivity,  is  given  by  the  following  equation, 

(Ref.  1): 
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•  • 


(11) 


where 


9(hK)=  f 


e  C  CV 


B 

K 

n 

e 

e 

c 

CO 

r 


is  the  magnetic  field  strength  . 

is  the  loading  factor  which  is  the  ratio  of  the  load 
voltage  to  the  open  circuit  voltage. 

is  the  electron  density. 

is  the  electron  charge. 

is  the  velocity  of  light. 

is  the  electron  cyclotron  frequency. 

is  the  average  time  between  electron-non-electron 
collisions . 

is  the  plasma  electrical  conductivity. 
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The  power  density  in  a  segmented  electrode  generator  is  given  by: 


■  ■  ■  (,2) 


where 

W  is  the  power  density. 

K  is  the  loading  factor  which  is  the  ratio  of  the  load 
voltage  to  the  open-circuit  voltage. 

is  the  plasma  electrical  conductivity. 

M|  is  the  gas  velocity. 

3  is  the  magnetic  field  strength  . 

The  calculations  are  performed  for  the  case  where  viscous  boundary  layer  forma¬ 
tions  are  negligible,  and  is  constant. 
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Sample  Results: 


MHD  generator  performances  are  presented  for  pure  lithium,  sodium, 
potassium,  rubidium,  and  cesium,  under  the  following  conditions: 

Total  pressure  of  10,  20,  and  30  psia 

Total  temperature  of  1400,  1600,  1800,  and  2000°K 

Av^\ach  number  of  0.8,  1.0,  1.2,  1.4,  1.6,  1.8  and  2.0 

Electron  temperature  range  of  static  gas  temperature  up  to  4000®K. 

Sample  results  are  shown  in  figures  1  -  13. 
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Discussion: 


Theoretical  performance  calculations  are  presented  for  lithium,  sodium, 
potassium,  rubidium,  and  cesium,  in  figures  1-5  respectively.  These  results  are 
based  upon  the  largest  elastic  electron-neutral  collision  cross  sections  reported  in 
the  literature.  In  the  case  of  rubidium,  a  linear  variation  between  potassium  and 
cesium  was  assumed.  A  comparison  of  alkali  metal  vapor  systems  is  presented  in 
figure  6.  A  complete  comparison  at  a  total  pressure  of  20  psia  and  a  total  temper¬ 
ature  of  1400°K  was  impossible  because  of  the  low  vapor  pressure  of  lithium . 

These  results  show  that  the  decrease  in  the  ionization  potential  of  the  heavier 
elements  is  greatly  overshadowed  by  the  increase  in  atomic  weight  and  especially 
by  the  increase  in  electron-neutral  collision  cross  section.  Lithium  is  undesirable 
because  of  the  high  operating  temperatures  required.  Rubidium  and  cesium  are 
undesirable  because  of  their  high  atomic  weights  and  large  electron  collision  cross 
sections.  Either  potassium  or  sodium  would  be  acceptable  working  fluids;  however, 
the  amount  of  equilibrium  dimer  In  the  vapor  is  much  greater  in  sodium  than  in 
potassium,  as  shown  in  figure  17.  Consequently,  potassium  appears  to  be  the  best 
choice  among  the  alkali  metals  as  an  MHD  working  fluid. 

Recent  data  (References  (8)  and  (9))  indicate  the  collision  cross  sections  for 
potassium  and  cesium  to  be  much  lower  than  has  usually  been  assumed.  If  this  is 
true,  the  MHD  generator  performance  is  greatly  enhanced,  as  shown  in  figures  7, 

8,  and  9.  The  influence  of  the  magnetic  field  upon  the  plasma  electrical  conduc¬ 
tivity  is  shown  in  figure  10.  The  influence  of  the  Mach  number  upon  the  performance 
is  shown  in  figure  1 1  .  The  effect  of  electron-potassium  collision  cross  section  is 
shown  in  figure  12.  The  relationship  between  the  plasma  electrical  conductivity 
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and  the  electron  temperature  is  shown  in  figure  13.  Thus,  if  for  potassium 
15  2  3 

Qen=  10  ■  cm^,  power  densities  of  1  kw/cm'^  can  be  expected  for  potassium  with 

a  total  gas  temperature  of  1400®K,  a  total  pressure  of  20  psia,  a  Mach  number  of  2, 
and  a  magnetic  field  around  64,000  gauss  with  a  loading  factor  of  0.75. 

For  convenience,  the  details  of  the  calculation  procedure  are  presented  in 
appendix  A.  Vapor  pressure  curves  of  the  alkali  metals  are  presented  in  appendix 
B.  The  vapor  pressure  data  presented  herein  were  obtained  from  R.  E.  Honig  (Ref,  3). 
Honig's  vapor  pressure  data  review,  the  most  comprehensive  collection  currently 
available,  is  based  primarily  upon  data  reported  in  References  (4),  (6),  and  (11). 
Curves  of  the  equilibrium  dimer  concentration  in  pure  alkali  metal  vapors  are  pre¬ 
sented  in  appendix  C.  As  expected,  the  relative  equilibrium  concentrations  of 
dimer  are  also  reflected  by  the  heats  of  dimer  dissociation  at  O^K,  Ref.  (2). 

Heat  of  Dimer  Dissociation 
Alkali  Metal  at  O^K 


*-'2 

1 .02 

e.v . 

Na2 

.73 

e.v . 

.514 

e.v. 

Rb2 

.49 

e  .V . 

Cso 

.45 

e.v. 

The  computer  routine,  used  to  calculate  the  MHD  generator  performance,  is 
presented  in  appendix  D. 
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APPENDIX  A 


Calculation  Procedure 
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Calculation  Procedure 


A.  The  quantities  to  be  specified  apriori  are: 

Tq,  the  total  temperature  of  gas, 

Pq,  the  total  pressure  of  gas,  atm 

My  the  Mach  number  in  the  MHD  generator,  dimensionless 

Tg,  the  electron  temperature  at  the  entrance  of  MHD  generator 
OK 

if ,  the  ratio  of  heat  capacities,  (cp/c^)/  dimensionless 

E  ,  the  first  ionization  potential  of  the  gas,  electron  volts 
o 

-  2 

Q  ,  the  electron-neutral  collision  cross  section,  cm 
en ' 

The  quantities  y,  and  are  determined  from  the  particular  gas  chosen. 

B.  Calculation  of  ,  (the  static  temperature  of  neutral  particles  in  the  MHD 
generator,  °K). 

The  quantity,  T, ,  is  determined  from  the  steady  isentropic  flow  equation: 

•  •  •  (13) 

where: 

T  is  in  OK 
o 

P,  is  in  atm. 

P  is  in  atm . 
o 

y  is  dimensionless 
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C.  Calculation  of  P|  ,  (the  static  pressure  of  gas  in  the  MHD  generator,  atm.). 


The  quantity,  Pi  ,  is  determined  from  the  steady  isentropic  flow 

equation; 


P. 


where: 


P 

o 


ar  - 1  M,  ^ 
2 


«7r  - 1 


atm . 


P  is  in  atm. 
o 

Y  is  dimensionless 
M,  is  dimensionless 


•  •  • 


(14) 


D.  Calculation  of  n^^,  (the  neutral  particle  density,  ^  of  neutrals/cm^ ) 


The  quantity  n^  at  the  entrance  of  the  MHD  generator  is  determined  from 

the  ideal  gas  law; 

n  =  P. 
n 


k  T. 


since 


k  =  1.38042  X  10  erg^^K,  (the  Boltzmann  Constant) 


1  erg  =  1  dyne-cm 

1  atm.  =  1 .013246  x  10^  dynes/cm^ 
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Thus 


n 


n 


:  f  1  atm . 

\  1 

r  1 

(  1.38042  x10"’^ 

1  dyne  -  cm  J 

= 


0.7340  X  1022  p, 


1  .013246  X  10  dynes 


cm- 


2  - 


atm . 


-] 


^  of  particles 

.3 


•  •  • 


(15) 


cm'- 


where: 

P|  is  in  atrr , 
is  in  OK 

E.  Calculation  of  ,  (the  electron  density,  ^  of  electrons/cm^) 

The  electron  density  is  determined  from  the  Saha  equation  based  upon 
the  electron  temperature: 


=  f2jrmg  g  ■  ^o/kTe 


Since 

”  9.1085  X  10”^®  grams  (electron  mass) 
k  =  ]  .38042  X  10”^^  ergs/^K  (the  Boltzmann  Constant) 
h  =  6.6252  X  10"^^  erg-sec  (the  Plank  Constant) 
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Thus 


0 


28 


3/2 


(2fl)(9.1085  X  10  gr.)  x 


1  .38042  X  10"’^  2^  ok] 


(6.6252  X  10  erg-sec.)^ 


3/2 


=  2.4146x10 


.+  15 


3/2  3 

erg  /  -  sec 


=  2.4146x10 


15  1 


cm 


Also  since 

1  e.v.  =  1  ,60207  x  10  joules 
1  joule  =  10.^  ergs 


Thus 


r 


(1  e.v.)  11.60207  X  10' 


-19 


( 1  .38042  X  10"’^ 

\  OK 


where 


=  11,606  Jo 
Te 


dimensionless 


E  is  in  e  .V, 
o 


T  is  in  OK 
e 
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Finally 


n  -  n 

n  e 


=  2.4146  X  10  e  ^o/^e  ,  ^  of  particles 

i3 


cm ' 


wfie 


(16) 


re 


is  ^  of  e!cctrons/cm^ 
is  ^  of  neutral^ cm^ 


Tg  is  in  OK 


is  in  electron  volts 
o 


F.  Calculation  of  ^  ,  (the  ratio  of  the  Debye  Shielding  length  to  the 
average  impact  parameter,  dimensionless) 

The  quantity  /\  is  determined  from  the  following  Equation  (Ref.  10): 

1/2 

/  l3  T  3  '' 

A  = _ 1 


2ZZ.e3 


f— 


Since  Z  =  Zj  =  1  in  this  cose 

e  =  4.80288  X  10  e.s.u,  (electron  charge) 

k  =  1.38042  X  lO”^*  eig/OK 

2  _ 


1  erg  =  1  gr  -  cm^/  sec^  =  1  dyne- 
1  e.s.u.  =  (dyne-cm^)  ^ 


cm 
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Thus 

A 


^1.38042  okJ 


(4.80288  X  10"’°e.s.u.)^ 


l/2l  3 


4 

=  1  .2388  X  10 


dimensionless 


(17) 


where 


T  is  in  OK 
e 

is  in  ^  of  electrons/cm^ 


G.  Calculation  of  QT  ,  ,  (the  conductivity  due  to  the  ions,  mhos/meter) 


The  quantity  OT  is  calculated  from  the  following  equation  (Ref.  10): 


tr  = 


-2  3/2 

1.51x10  T^ 


^n  A 


mho 

meter 


(18) 


where 


Tg  is  in  OK 


A  's  dimensionless. 


H.  Calculation  of  ^ ,  (the  conductivity  due  to  the  neutrals; 
mhos/meter) 

The  quantity  is  determined  from  the  scalar  conductivity  equation 

(Ref.  1); 

n  ■ 

m  0 

where 


mg  =  9.1085  X  10  grams 

k  =  1 .38042x1  o'’ ^ergs/®K 

1  2/2 
1  erg  =  gr.  crn  /  sec 

,  _  1/2  3/2/ 

1  e  .s.u .  =  gr  cm  /  sec 

cm /sec  =  (1/9x10"  ohm) 
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Thus 


-10 

(4.80288  x  10  e.s 


/% 

V- 


n 


-  "e)  Qen 


•) 


-28 

(9.1085  x  10  gr)  /  J_ 


1  .38042  ^  °K 
OK 


,1/2 


TT  9.1085  X  lo'^®  gr 


cm  V2/„; 


e .s.u . 


-12  r 

c 

^  L 


1 


9  X  10 


1 1  ohm-cm 


1  r  J0^_cm_l 
J  I  mefer  J 


sec 


^  4.525  X  10'®  n 

0^  =  _ L 


mho 


T  1/2,  V  ^  meter 


^  *  (19) 


e'  en 


wh 


ere 


T  is  in 

e 

is  ^  of  neutrals/ cm  ^ 
n  is  ^  of  electrons/ cm^ 


Qen  is  in  cm 


I.  Calculation  of  (the  plasma  electrical  conductivity,  mhos/  meter) 


The  plasma  electrical  conductivity  Is  determined  from  the  following 


equation: 


-1 


-1 


r,  "rn 


or 


07  (t: 


irr  *(77 


mhos 

meter 


(20) 
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yyhere 


IS  in  mhos/ meter 
^  is  in  mhos/ meter 

J  .  Calculation  of  (dimensionless) 

The  quantity  6/7^  necessary  to  sustain  the  given  is  determined  from 

the  equation  as  shown  below  (Ref.  5): 

1/2 

^  -1 
T. 


coT  = 


dimensionless 


where 


T  is  in 
e 

T,  is  In 

^  is  dimensionless 
Ml  is  dimensionless 

K  .  Calculation  of  B(1-K)  (the  magnetic  field,  gauss) 

The  quantity  B(l-K)  (the  magnetic  field  necessary  to  sustain  the  desired 
when  the  generator  is  short  circuited),  is  calculated  from  the  equation 


B(NK)  = 


r  %ec Cor  \ 

\  <^  ) 


The  obove  equation  is  derived  from  the  following  simple  relations  (Ref.  1): 
CO  =eB(1-K)/c 
and 
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Since 


e 

c 

1  e  .s .u . 
1  gauss 
1  sec 


=  4.80288  X  10  e,s.u.  (electron  charge) 
=  2  .998  X  10^^  cm/ sec  (speed  of  light) 

=  9,'^2/cm'/2-»c 
=  9x10"  ohm-cm 


Thus 


B(NK) 


of  electrons\/4. 80288  x  10  ^^e.s.u.  j  / 2  993  x  10^^-^)  ^ 

electron  /  V  secy 

mho 

meter 


r  1/2  3/2  Tr- 

'  gr  cm  1 1  gauss 


sec 


e.s.u . 


sec 


[1  meter  T  f  l/ohm  T  f* 


1 1  ohm-cm 


sec 


J 


,  gauss 


Thus 


B(l-K)  =  1.600  X  10' 


.9  ^  £vr 


(22) 


wh 


ere 


K 


ct,r 


is  dimensionless 
is  ^  of  electrons/cm^ 
is  dimensionless 
is  in  mhos/meter 
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L,  Calculation  of  U,,  (the  gas  velocity  at  4ie  entrance  of  the  MHD 


I 

5 

generator,  meters/ second) 

The  quantity  is  determined  from  the  sonic  velocity  equation: 


Since 

R  =  1 .987  cal/gr.  mole  (the  gas  constant) 

1  calorie  =  4.186  joules 

1  joule  =  lO^ergs 

2  2 

1  erg  =  1  gr-cm  /  sec^ 

Thus 
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where 


T. 


is  dimensionless 
is  in 


is  at,  wt.,  gr./ gr.  atom 
M,  is  dimensionless 

M.  Calculation  of  (the  MHD  generator  power  density,  kilowatts/cm'^ ) 

'y  "  ■ 

Thus  2 

W  =  (T  (1-K)  K  U,B 
P 


or 


W(l-K)  = 
K 


^  U.^Bd-Kj 


Since 


Thus 


10  gauss  =  1  weber/ meter  ^ 
1  weber  =  1  volt  second 


W(1-K)  =  (T^U.  [b(1-K)]' 

K 


r  weber/ meter^ 
10^  gauss 

1  meter  ^ 


2  2 

mho-meter  gauss 

meter  -  sec  ^ 

watt-ohm 


f 


lo'*  cm^ 


^*7  r  volt-sec .  *7  r  1  /ohm  1  T wott-ohm  *1 

•J  I,  weber  JL  mho  J  L  volt^  J 
JL  10^  watts  J 


W{l-K)=  10 
K 


watts 


KW 

cm 


(24) 
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where 


K  is  dimensionless 

(TV  t.in  mhos/meter 

Ui  is  in  meters/second 

B  is  in  gauss 


APPENDIX  B 


Vapor  Pressure  Curves  for  Alkali  Metals 
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Vapor  Pressure  of  Alkali  Metals 


The  vapor  pressure  data  presented  herein  were  obtained  from 


R  E.  Honig  (Ref.  3),  Vapor  pressure  curves  are  shown  in  figures  14-15. 
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APPENDIX  C 


Equilibrium  Dimer  Concentration  in  Pure 


Alkali  Metal  Vapors 
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Equilibrium  Dimer  Concentration  in  Saturated  Pure  Alkoli  Metal  Vapors: 


Consider  the  following  reaction: 


2M 


(25) 


where  M  is  an  alkali  metal .  As  usual,  the  equilibrium  constant  for  the 


homogeneous  chemical  reaction  is  defined  as 

K  =(Pm)^ 

P 


(P  M2) 


.  .  •  (26) 


where  is  the  partial  pressure  of  monomer  in  the  vapor 
is  the  partial  pressure  of  dimer  in  the  vapor 


The  mole  fraction  of  dimer  is  given  by  the  ratio 


X2  =  P 


M2 


•  •  • 


(27) 


where  P^  is  the  static  pressure  of  the  vapor 
Consequently, 


X2  =  (2Po  +  K) 


+  4K  P 
P  P  o 


2  P. 


The  equilibrium  constant  is  obtained  from  the  well-known  equation 


(28) 


(29) 
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where 


is  the  absolute  free  energy 

R  is  the  universal  gas  constant 

T  is  the  static  temperature 

refers  to  the  monomer 
X2  refers  to  the  dimer 

Thus  a  knowledge  of  the  static  pressure,  static  temperature  and  free  energy  permits 
the  dimer  concentration  to  be  calculated  from  equation  (28). 

Such  calculations  have  been  performed  by  AAeisI  and  Shapiro  (Ref.  7). 

It  should  be  noted  that  the  free  energies  of  the  alkali  metals  were  obtained  from 
Reference  1  .  Dimer  curves  are  shown  in  figures  16-18. 
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APPENDIX  D 

Program  Flow  Diagram 
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PROGRAM  FLOW  DIAGRAM* 


♦ 

The  numbers  in  parentheses  refer  to  the  equations  as  numbered  in  Appendix  A. 
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APPENDIX  E 
Physical  Properties 
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Physical  Properties  Used  for  MHD 


Performance  Calculations 


(1) 


Element 

0) 

•  m 

<  o 

2 

em 

m 

n 

gr/gr.  atom 

Lithium 

5/3 

5.363 

2  X  lO'^"^ 

6.94 

Sodium 

5/3 

5.12 

3  X 

22.997 

Potassium 

j 

5/3 

4.318 

4x  10"’"^ 

3  X  10"J5 

1  X  10"’^ 

39.1 

Rubidium 

5/3 

4.159 

4.7  X  lO"’"^ 

85.48 

Cesium 

5/3 

3.87 

5.3  X  ]0~\^ 
3.6  X  10"'^ 

132.91 

(1)  Calculations  were  performed  for  old  literature  values  of  and  also 
for  the  latest  literature  values,  (see  references  4  and  5). 
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APPENDIX  F 
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Table  of  Nomenclature 


T 

e 

To 

y 

CO 

r 


'”1 

€ 

K 


o 

rr 


A 


(Listed  in  the  Order  of  Appearance) 

Electron  Temperature, 

Total  gas  temperature, 

C 

Ratio  of  heat  capacities,  (  p  ),  dimensionless 
Electron  cyclotron  frequency,  1  /  sec. 

Average  time  between  electron-nonelectron  collisions,  sec. 

Mach  number  in  MHD  generator,  dimensionless 

Loss  Factor,  to  account  for  inelastic  collisions,  dimensionless 

Loading  Factor,  ratio  of  load  voltage  to  open  circuit  voltage, 
dimensionless 

Electron  density,  ^  of  electrons  /  cm^ 

u  3 

Neutral  particle  density,  "  of  neutrals  /  cm'^ 

First  ionization  potential  of  vapor,  e.v. 

Static  pressure  of  gas  in  MHD  generoror,  atm. 

Static  temperature  of  neutral  particles  in  MHD  generator, 

Total  gas  pressure,  atm. 

Total  gas  temperature, 

Electrical  conductivity  due  to  ions,  mhos/ meter 

Ratio  of  Debye  shielding  length  to  the  average  impact  parameter, 
dimensionless 

Electrical  conductivity  due  to  neutrals,  mhos/ meter 
Scalar  electrical  conductivity  of  plasma,  mhos/ meter 
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u, 


'"r 

W 

Q. 

V 

_E_ 

M 

K 


en 


M„ 


F° 

k 


m  - 


h  = 


e  = 


R  = 


Velocity  of  plasma  in  MHD  generator,  meters/ second 
Atomic  weight  of  neutrals,  gr./ gr.  atom 
Power  density  in  MHD  generator,  kw/cm'^ 

2 

Electron-neutral  collision  cross  section,  cm 

Mean  electron  speed,  (Maxwellian  distribution  assumed),  cm/ sec. 

Electric  field  vector  =  i  E  +iE  +kE 
—  X  y  —  z 

Current  vector  =  j  J  y  j  J  y  1^  Jz 
Any  member  of  the  alkali  metal  series 
Equilibrium  constant  defined  by  equation  26,  atm 
Partial  pressure  ofmc  iomer,  atm 
Partial  pressure  of  dimer,  atm 

Absolute  free  energy,  cal./gr,  mole 

1  .38042  X  10  ergs/oK  (Boltzmann's  constant) 

-28 

9,1085  x10  grams  (electron  mass) 

6.6252  X  10"^^  erg-sec,  (Plank's  constant) 

4.80288  X  10  ^^e.s.u.  (electron  charge) 

2.998  x10^  cm /sec.  (speed  of  light) 

1 .987  cal./gr.  mole  -®K  (gas  constant) 
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Conversion  Units 


1  -tm  = 

1.013246  X  10*^  dynes/ 

1  erg  = 

1  dyne  -  cm 

1  e.v*  = 

1.60207  x10"’’  joules 

1  joule  = 

10^  ergs 

1  e.s.u .  = 

2.  1/2 

^dyne  -  cm  ; 

1  dyne  = 

/  2 

gr.-cm/sec 

1  ohm  = 

/  1  .  sec/cm 

‘  9xl0"> 

1  gauss  = 

gr.  '/V  c.  1/2.  „c. 

1  calorie  = 

4.186  joules 

2 

1  weber/meter  = 

1 0  gauss 

1  weber  = 

1  volt-sec . 
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Figure  6.  Comparison  of  Power  Densities  for  Pure  Alkali  Metal  Vapor  Systems 
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Figure  7.  Power  Density  for  Cesium,  Q  =  3.  6  x  10 


B(l-K), GAUSS 


potassium:  y=5/3 

Eo=  4.318  t.v. 
Q,n=3xlO-'»cm 
Pq  =  20.0  p.t.i.a. 

- To=  1400*  K 

- To  =  1600*  K 

- To=  1800*  K 

- To=  2000*  K 


Figure  8, 


Power  Density  for  Potassium,  Q  =3x10 


Figure  11*  Influence  of  Mach  Number  for  Power  Density  of  MHD  Generator 
for  Potassium  Vapor 


Figure  13.  Relation  Between  Plasma  Conductivity  and  Electron  Temperature  for  Potassium  Vapor 


Figure  15.  Vapor  Pressure  Curves  of  Alkali  Metals 


Figure  16.  Equilibrium  Dimer  Concentration  in  Saturated  Pure  Alkali  Metal  Vapors  [  Calculated  from  the 
Equilibrium  Molecular  Weights  as  Reported  by  Meisl  and  Shapiro,  ref.  (7)  ] 


Figure  17.  Equilibrium  Dimer  Concentration  in  Pure  Alkali  Metal  Vapors 
Near  Atmospheric  Pressure 
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Figure  18.  Equilibrium  Dimer  Concentration  in  Pure  Potassium  Vapor 
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